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Hybrid graphene-topological insulator (TI) devices were fabricated using a mechanical transfer 
method and stndied via electronic transport. Devices consisting of bilayer graphene (BLG) under 
the TI Bi 2 Se 3 exhibit differential conductance characteristics which appear to be dominated by 
tunneling, roughly reproducing the Bi 2 Se 3 density of states. Similar results were obtained for 
BLG on top of Bi 2 Se 3 , with 10-fold greater conductance consistent with a larger contact area 
due to better surface conformity. The devices further show evidence of inelastic phonon-assisted 
tunneling processes involving both Bi 2 Se 3 and graphene phonons. These processes favor phonons 
which compensate for momentum mismatch between the TI T and graphene K, K' points. Finally, 
the utility of these tunnel junctions is demonstrated on a density-tunable BLG device, where the 
charge-neutrality point is traced along the energy-density trajectory. This trajectory is used as a 
measure of the ground-state density of states. 

PACS numbers: 73.40.Gk, 73.50.-h, 73.20.At 


In topological insulators (TIs), the protected surface 
state is bound to the interface between materials with dif¬ 
ferent bulk topological invariants [T]. Although the sur¬ 
face state will remain gapless as long as the symmetries 
protecting the bulk topological invariant are retained, its 
dispersion and location depend on the band-structures 
of the two interfacing materials. A model system allow¬ 
ing the investigation of such interface properties is the 
graphene-TI hybrid, where graphene resides in immedi¬ 
ate proximity to the TI, separating it from the vacuum 
or from another trivial dielectric. Indeed, as part of a 3D 
continuum, graphene is topologically trivial. However, 
with a significant spin-orbit (SO) term, graphene could 
become a 2D TI [5]. 

Recent studies nn which theoretically investigate 
the properties of such hybrids, generally assume that 
the graphene and TI bands couple strongly, leading to 
fundamental modifications to the graphene band struc¬ 
ture, which is expected to inherit an enhanced SO cou¬ 
pling j3] and attain non-trivial spin-textures [6]. When 
sandwiched between two ultra thin layers of Bi 2 Se 3 , the 
graphene layer may become a 2D TI [^. However, as the 
surface states of TIs such as Bi 2 Se 3 are centered at the T 
point in fc-space, they should be prevented by momentum 
conservation from hybridizing with the graphene bands, 
centered at the K, K' points. This difficulty is often re¬ 
solved theoretically by assuming that periodicity at the 
interface would cause the graphene Brillouine zone (BZ) 
to fold (SM]: bringing significant spectral weight to the 
r point. 

Here we present an experimental study of graphene- 
TI hybrid devices fabricated using a mechanical transfer 
technique [7] which allows the fabrication of vertical het¬ 
erostructures involving a wide variety of van-der-Waals 
materials. In our devices, which consist of either mono- 



FIG. 1: (a) Optical image of Device 1, consisting of BLG on 
the bottom and BbSes on top. BLG outline marked in green, 
(b) Device 1 with contacts. Bar is 2/rm long, (c) dl/dV vs. 
Vti of Device 1 (source on the Bi 2 Se 3 .). Inset: Schematic 
showing the BLG (green) underneath the TI (teal, outline 
representing the TI surface states), (d) Device 2: Single BLG 
flake deposited on top of two Bi 2 Se 3 flakes (e) AFM phase 
image of Device 2, arrows mark the two separate junctions (f) 
dl/dVvs. Vs_D of the two junctions in Device 2. Top junction: 
blue, bottom: red. Inset: Schematic showing graphene or 
BLG on top of the TI. (g,h) Annotations illustrating the DOS 
alignment in tunneling between BLG (g) and MLG (h) and a 
TI. All dl/dV data are normalized per fim^. 

layer (MLG) or bilayer graphene (BLG) stacked above 
or below the TI Bi 2 Se 3 , the graphene and TI layers are 
weakly coupled, with transport between them governed 
by tunneling. The devices reveal a rich inelastic spectrum 
consisting of graphene and Bi 2 Se 3 phonons. Finally, we 
demonstrate the utility of the interface tunnel-junction 
to probe the density-dependent evolution of spectral fea¬ 
tures in BLG. 
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Fig. 1(a) shows Device 1, where a Bi 2 Se 3 flake is 
transferred on top of a BLG flake (outlined in green) de¬ 
posited on Si02. The flakes are independently contacted 
by evaporated metallic electrodes (Fig. 1(b)). The BLG- 
TI junctions are studied by measuring the differential 
conductance dl/dV vs. applied bias voltage Vsd with 
the source on one material and the drain on the other 
(Til to Gl). A four-probe geometry is realized by mea¬ 
suring the voltage across the opposing contacts {Vti = 
V(TI2)-V(G2)). The dl/dV trace at T = AK of Device 

1 (Fig. 1(c)) shows a pronounced suppression at nega¬ 
tive Vti values, similar to STM measurements taken on 
Bi 2 Se 3 [8HT0]. where dl/dV is proportional to the density 
of states (DOS) and the suppression at negative bias is 
associated with the bulk gap of the Bi 2 Se 3 . The BLG-TI 
interface therefore behaves as an effective tunnel junction 
although no intentional barrier was placed between the 
materials. It is a remarkably robust junction, maintain¬ 
ing a stable signal while sustaining high voltage biases 
exceeding 0.5 V at negative bias. 

The tunneling functionality of the graphene-Bi 2 Se 3 in¬ 
terface could be a consequence of Bi 2 Se 3 oxidation m- 
In over 20 devices studied, the interface resistance varied 
from 10 to 10 kil^rrV, which could be associ¬ 

ated with varying degrees of oxidation. The stability and 
high bias achieved by graphene-TI junctions, however, do 
not favor this explanation - oxide-based junctions rarely 
function at biases in excess of 200 meV [12]. Alterna¬ 
tively, the variation in interface resistance could arise due 
to differences in the effective contact area: graphene con¬ 
forms to the underlying substrate and has angstrom-scale 
height variations when deposited on Si 02 resulting in an 
effective partial contact area. To test this we studied de¬ 
vices where BLG is deposited on top of Bi 2 Se 3 : in Device 

2 a single BLG flake covers two Bi 2 Se 3 flakes (Fig. l(d- 
e)). Measured separately (red and blue curves in panel 
(f)), the two junctions exhibit comparable dl/dV char¬ 
acteristics which are very similar to the “graphene on 
bottom” devices (albeit with smaller bias range). 

The interface conductivity of “graphene-on-top” junc¬ 
tions is an order of magnitude greater than “graphene on 
bottom” ones, indicating that the effective contact area 
is a plausible source of variation. It is not clear, however, 
what is the origin of the generic tunneling functionality. 
An interesting possibility is that the tunnel barrier is as¬ 
sociated with the lack of direct chemical bonding between 
the two layers. Inter-layer tunnel barriers in layered ma¬ 
terials is observed in high Tc superconductors [13] and 
in vdW materials such as 4Hb-TaS'2 M- Incoherent in¬ 
terlayer transport was also reported in stacked twisted 
bilayer graphene devices |15j . and graphene layers could 
be weakly coupled to underlying graphite m- Testing 
this hypothesis would require further investigation. 

As expected, the dl/dV characteristics vary with the 
choice of MLG vs. BLG. The tunneling current is ex¬ 
pressed as a convolution of the graphene DOS pai^) and 
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FIG. 2: (a) dl/dV vs. Vti of Device 3: BLSea on top of 
monolayer graphene. Inset: Closeup of the data showing the 
zero-bias suppression associated with inelastic tunneling fea¬ 
tures. (b) d?I/dV^ vs. Vti of Device 1 (BLG) and Device 2 
(MLG). Both devices exhibit a peak in the 2nd derivative, cor¬ 
responding to an inelastic excitation, at 7.7 meV, associated 
with a BbSes optical phonon (c) dl/dV vs. Vti of Device 4, 
which is a vertical heterostructure of a hBN substrate, BLG, 
and BbSea on top. Inset: Image of the device (Bar is 4jj.m). 
(d) Id^I/dV^l vs. Vti of Device 4, exhibiting a rich spec¬ 
trum of inelastic features. The energies marked in the figure 
are discussed in the text. Inset: Phonon assisted tunneling 
in momentum mismatched materials: Tunneling from the TI 
states at the F point to graphene states at the K, K' points 
is assisted by BZ-edge phonons. 


the TI DOS PTii^) [T7] 


I 



eV)pTi{s)\t{e)\'^de 


( 1 ) 


where A is the effective overlap area and t{e) stands for 
the tunnel coupling through the barrier. This expression 
is qualitatively depicted in the schematics in Fig. l(g,h), 
showing pG on the left and pri (including bulk and sur¬ 
face) on the right. For BLG devices (panel (g)) pai^) 
is relatively featureless and can be factored out of the 
integral, so the dl/dV curve effectively probes pTii^)- 
This explains why BLG-TI dl/dV traces are similar to 
STM measurements of Bi 2 Se 3 , and suggests that BLG 
may be useful as a tunneling electrode for probing other 
vdW systems. MLG-TI devices (panel (h)) have energy- 
dependent p(e) on both sides of the barrier, and Device 
3 (Fig 2(a)) indeed exhibits a very different dl/dV char¬ 
acteristic, with a stronger suppression around zero-bias 
which we associate with the graphene Dirac point. 

Tunnel junctions are useful in measuring inelastic spec¬ 
tra, where the onset of processes such as phonon-assisted 
tunneling appears as step-increase features in the differ- 
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FIG. 3: (a) Schematics of gate-dependent tunneling and transport processes. (I) BLG at n = 0, Vti at zero bias. (II) tracing 
the CNP in a transport measurement: Vbg and Vtg are balanced to keep n = 0. (Ill) tunneling into the GNP: Vbg < 0 
tunes the BLG to n > 0, Vti < 0 keeps tunneling to the CNP (IV) phonon-assisted tunneling into the CNP. (b) d^I/dV^ vs. 
Vti obtained from integrating data in (e). (c) In-plane electronic transport vs. Vbg and Vtj. Annotation II marks the CNP 
diagonal trajectory. Annotation I marks the point where this trajectory crosses Vti = 0. (d) dl/dV vs. Vti at a range of gate 
voltage (5V interval), (e) d?I/dV^ vs. Vti and Vbg- Annotations mark the elastic process (III) and phonon-assisted processes 
(IV). The lines superimposed on the data are fits to the trajectories of elastic and inelastic tunneling into the CNP |18| . 


ential conductance. These are detectable as peaks in the 
2nd derivative d^I/dV^ at bias voltages corresponding 
to the phonon energies. It is possible to differentiate be¬ 
tween Bi 2 Se 3 and graphene phonons because their respec¬ 
tive spectra span different energy ranges (up to 20 meV 
for Bi 2 Se 3 [19], and up to 200 meV for graphene [20]). 
Several devices exhibit well-developed inelastic spectra: 
Device 3 (Fig. 2(a)), has sharp step-like features in 
dl/dV close to zero (inset). The corresponding dtd’IjdV'^ 
plot (Fig. 2(b)), shows peaks at ±7.7 meV which coin¬ 
cide in Device 3 (MLG) and in Device 1 (BLG). We as¬ 
sociate this feature with a Bi 2 Se 3 surface optical phonon 
previously identified using Helium scattering at the same 
energy |21] . This phonon plays a role in suppressing sur¬ 
face transport |22j . and inelastic tunneling data may be 
useful in probing its coupling to surface electrons. 

Signatures of graphene phonons, appearing at higher 
bias voltages, are found in numerous devices. For exam¬ 
ple, in Device 4 (Fig. 2(c,d)), a higher quality device 
fabricated by deposition of BLG on h-BN followed by 
transferring Bi 2 Se 3 on top. The d?I/dV"^ plot shows a 
prominent peak at 67 meV, corresponding to the energy 
of the graphene ZA/ZO mode at the K,K' points. This 
phonon is seen in many graphene and graphite tunneling 
experiments both in devices [23] and in STM PHHS], and 
is generally believed to assist the tunneling process by 
providing the momentum required to inject an electron 


to the K, K' points. In our devices, where graphene and 
Bi 2 Se 3 Fermi surface momenta are highly mismatched, 
it is likely that the same phonon is required to allow 
momentum-conserving tunneling from low energy Bi 2 Se 3 
states, centered at the F point, to the graphene K,K' 
points (Schematic in inset to Fig. 2(d)). Device 4 also 
exhibits a feature at 162 meV associated with the LA/LO 
mode, also a K, K' phonon, and a feature at 195 meV as¬ 
sociated with the TO mode at F. Other features (e.g. a 
phonon at 36 meV) could be associated with the h-BN 
substrate [23] . 

In the rest of this Letter we focus on Device 5, which 
has the same geometry as Device I: Tl on top of a 
density-tunable BLG flake. Device 5 can be measured 
by graphene in-plane transport (G1 to G2) or by tun¬ 
neling (Tl to GI). In-plane resistance R 2 t{Vti,Vbg) 
is presented in Fig. 3(c), showing a resistance peak 
when the chemical potential crosses the BLG CNP. The 
Bi 2 Se 3 electrode acts as a well-behaved gate, and the 
small tunneling current does not interfere with the in¬ 
plane measurement. The resistance map is typical of 
doubly-gated graphene [Ml HZ], where the diagonal fea¬ 
ture corresponds to a high resistance state in the BLG re¬ 
gion underneath the TL This feature intercepts Vti = 0 
at Vbg = Fq = 20V, indicating that the BLG is p-doped 
and the CNP is energetically mismatched from the Tl 
DP (Fig. 3(a) annotation I). Its slope is dictated by the 
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requirement that the top and bottom gates compensate 
each other’s charge (annotation II) and therefore follows 
the ratio Cbg/Cti {Cti and Cbg are the graphene-TI 
and graphene-back gate capacitances, respectively). Us¬ 
ing this relation, we extract the geometric capacitance 
of the TI-BLG junction Cti = 1-3 x 10“^U/m^, 110 
times greater than the back-gate capacitance. Both the 
relative permittivity of the interlayer medium, e, and 
the graphene-Bi 2 Se 3 effective distance, d, are unknown. 
However, the measured value of Cti fixes their ratio to 
1.5 nm“^. As the interface is a tunnel junction, d is 
unlikely to exceed 3 nm, setting a limit of e < 5. 

Gate-dependent tunneling measurements are presented 
in Fig. 3(d-e): Such measurements are sensitive to energy 
shifts in the spectral features due to the changes in den¬ 
sity. In Fig. 3(d) it is clear that the underlying structure 
of the Bi 2 Se 3 gap at negative Vri is retained while a set 
of features closer to zero bias evolve with Vbg- To trace 
these features we plot the second derivative d'^I/dV'^ as 
a color map in Fig 3(e). The data exhibit a set of diag¬ 
onal features with opposite slopes and gate-independent 
features which appear as faint vertical lines. The lat¬ 
ter gate-independent features are enhanced by averaging 
over all back-gate values (Fig. 3(b)), where we again 
find the ZA/ZO phonon features at ±65 meV. The 22 
meV feature coincides with the energy of the Bi 2 Se 3 Afg 
phonon [TU] . 

Vbg changes the graphene density, vertically shifting 
the band-structure (Fig. 3(a), annotation III). To trace 
a spectral feature such as the CNF, a voltage SVri has to 
be applied to compensate for the density-induced change 
in chemical potential, dfiG ■ This is formulated as: 

cSVti = SfiG (2) 

To trace spectral features on the Vti-Vbg plane, we 
have to note that planar tunneling electrodes are large 
capacitors which charge the graphene layer at finite bias, 
as discussed earlier. The equations governing the charg¬ 
ing of the graphene layer connect the incremental gate 
voltages 5Vti,bg to the incremental charges on the back- 
gate and TI 6nBG,Ti' 

- eSn, = C, (^Vi - (3) 

where i = BC,TI. 

Interestingly the GNP tunneling condition in Eq. 
also ensures that there will be no extra charge accumu¬ 
lated on the TI {Stiti = 0), as the argument in paren¬ 
thesis in Eq. [^remains zero. More generally, this means 
that constant energy features lie on trajectories which 
keep the charge on the tunneling electrode hxed. In 
the Supplementary Information |I8] we derive the CNP- 
tunneling trajectory and find its slope to be: 

SVti _ Cbg 
SVbg Cq + Cbg 


where Cq = e^p is the quantum capacitance of the BLG. 
The slope is independent of Ctb consistent with the in¬ 
duced charge on the TI remaining zero. 

As a result, Eq. |^can be used to evaluate Cq by fitting 
the slope of the features appearing in the figure. The 
most prominent feature in Fig. 3(e), marked by “III”, is 
likely the charge neutrality point, which attains a cusp 
in the DOS due to a finite displacement field (annotation 
I). The slope of this feature is not constant, and appears 
to vary along the trajectory (Fig. SI). This indicates 
that Cq, and hence the DOS, changes with density, as 
expected for the hyperbolic dispersion of BLG [^SHSO]- 

We fit the expected trajectory (Fig. 3(e), feature 
“III”). To calculate this fit it is crucial to acknowledge 
that the BLG dispersion varies within the Vbg — Vti 
plane, where at each point the displacement field is dif¬ 
ferent. The calculation, described in the SI [TH], is carried 
out self-consistently, using an approximate trajectory for 
calculating the displacement held at each value of Vbg 
and Vti- We then extract n{e), p{n) and hnally Cq and 
integrate the trajectory from the slope in Eq. The 
calculation is rehned iteratively. Using the Eermi veloc¬ 
ity vp as a htting parameter, the full trajectory can be 
reproduced with vp = 1-06 x 10®to/s, in good agreement 
with non-interacting values for BLG. 

Multiple replica (marked by IV in Fig. 3(e)) of the 
GNP tunneling feature, appear as lines running parallel 
to the elastic feature (“IV”). They are associated with 
phonon-assisted inelastic tunneling processes to the GNP. 
Unlike the gate-independent phonon-onset features dis¬ 
cussed above, where tunneling takes place to the Fermi 
energy, these features represent phonon-assisted tunnel¬ 
ing to some sharp spectral feature (here the GNP), and 
evolve with the gate in parallel to the elastic tunneling 
feature. Their trajectories, presented in Fig. SI, should 
depend on the elastic feature, after accounting for mod¬ 
ified displacement helds at higher bias. However, to fit 
the actual data we find that the Fermi velocities have 
to be modified in each of the inelastic features (assum¬ 
ing interlayer coupling t± remains fixed). The modified 
Fermi velocities are 0.98 x 10® m/s for the -22 mV fea¬ 
ture, and 0.9 x 10® m/s for the -65 mV feature. This 
points, perhaps, to velocity renormalization. 

In summary, the graphene-Bi 2 Se 3 interface is a high 
quality tunnel junction which can be integrated in to a 
density-tunable device. Further studies are required to 
address the effect of junction properties such as crystal¬ 
lographic orientation, interface quality and the effect of 
Bi 2 Se 3 oxidation. Phonon-assisted tunneling is observed 
at finite bias, with specihc phonon-activation processes 
which bridge the mismatch in crystal-momentum be¬ 
tween the two materials. Nevertheless, tunneling in gen¬ 
eral is not momentum conserving, probably due to junc¬ 
tion inhomogeneity or other scattering processes. The 
role of the surface state in the tunneling process is also 
an open question: where some device (Device 1) seem 
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to reveal evidence of bulk states in the tunneling signal, 
others (Device 4) appear to be mostly surface-dominated. 
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